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ABSTRACT 


A general approach has been developed for computing 
the force due to solar radiation on an object of arbitrary shape. 
This method circumvents many of the existing difficulties in 
computational logic presently encountered in the direct analytical 
or numerical evaluation of the appropriate surface integral. It may 
be applied to complex spacecraft structures for computing the total 
force arising from either specular or diffuse reflection or even 
from non-Lambertian reflection and re-radiation. 
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SECTION 1 - INTRODUCTION 


The problem of computing the total force or total torque 
on a spacecraft due to solar radiation is, in general, veiy difficult. 
Mathematically, it requires the evaluation of a surface integral 
over only the illuminated region of the surface. Even if the illu- 
minated region is known by some other means, the evaluation of 
the surface integral can still be very difficult analytically 
in the case of complex spacecraft structures. Moreover, if the 
illuminated region is not known a priori, the difficulties are 
compounded by having to determine self-shadowing. For non-convex 
objects, it is not trivially governed by a condition such as 
cos 0 ^ 0 where Q is the angle between the sun vector and the out— 
wan'd vector normal to the surface. In fact, the logic in the present 
methods becomes extremely complicated and is also not fool-proof. 
Additional difficulties are introduced by choosing a set of points 
(vertices) on the surface to form a network in approximating it; 
this inadvertently leads to book-keeping problems associated with 
selecting appropriate sets of points for computing surface elements. 

This paper presents a general method for performing the 
computations without encountering the difficulties described above. 

It does not attempt to evaluate the surface integral directly as 
it presents itself as done in the usual methods, but considers 
the same problem from a slightly different point of view which 
leads to the same results. 
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For convenience, let us use the following notation; 
u “ unit vector along a specified direction 
V * any unit vector orthogonal to u, i.e., u • v ■ 0 
w * third unit vector forming orthogonal triad, i.e., w = u x 'v 
■ 0 = origin of coordinate system 
P “ any point on object's surface 
"r “ vector from 0 to P 
P'= projection of P onto (v,w)-plane 
projection of "r onto (v,w)-plane 
(x,y,z) * reference orthogonal system for describing object's 
surface. 

In the present anailysis, it is advantageous to choose u to be opposite 
in direction to the incident solar radiation. (Alternatively, it can 
also be chosen to be in the same direction.) 

The vectors u and "r are known in the (x,y,z) system. 

••A 

In general, if V is any vector, then it may be more explicitly 
written in the (x,y,z)-space as ^ and has components V^, V^, V^. 
That is, we implicitly mean 




0^0 


In view of the definition of the vector v, we may choose 
0. Then, it may be shown that the other two components are given by 
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Then, using equation (1.6), the vector "t/ ^ is transformed 

vx,y, zy 

and we obtain 


= 


-/• 


6 /. 7 ) 

/I,. = 





* 

An, 

V- 

A, «y -f AjKTj 

6A?) 


Consequently, the projection vector r* is simply given by 

a' m & -h Cj (I'IO) 

The component r^ of the vector "r is particularly important because, 
for a complex spacecraft structure, it can be used to yield the sxirface 
element which is directly exposed to solar radiation^ This can be seen 
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as follows: For any given point on the (v,w)-plane (i.e,, for any 

given vector'?'), the point on the spacecraft which is not shadowed 

is the one which has the maximum value of r^, independent of where 

the origin of the (u,v,w) coordinate system is chosen* (It would be 

the minimum value of r if the vector u had been chosen to be in 

u 

the same direction as the incident solar radiation.) To find the 
illuminated svirface of the spacecraft, we proceed by dividing 
the (v,w)-plane into cells of area Av Aw with cell centers (v,,w,). 

X J 

At these cell centers, the illuminated surface element is the one 
which has the maximum value of r^. In this way, the logic of deter- 
mining self-shadowing is extremely simple as compared to other methods 
which encounter considerable difficulty conceptually and computationally. 
Thus, given a vector ?' = (0,v^,w^), the vector? “ ^ ^ 

corresponding to the illuminated point is determined. It is then 
transformed to the (x,y,z)-space by the equation 





A, , > 


( A ^0 


At this point ?, the unit vector n, ■. normal to the surface 

(x,y,z)» (x,y,z) 


is then obtained by 




W\ 


0-n) 


where ^(x,y,z) = 0 denotes the equation of the surface in a region 
containing ?. For convenience, the direction of n is chosen such that 





This choice of direction automaticaQ-ly makes n the outward unit normal 
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if the surface element belongs to a closed surface. Moreover, 
it establishes a direction for n in the case of a surface for which 
an outward unit normal is meaningless (such as a finite planar surface). 
The vector n^^ ^ is then transformed to the (u,v,w)-space using 
the equation 


A 

■n 
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The cell (v. ,w. ) whose area is corresponds to 

a surface element whose area is denoted by AA. It is evident that 
we have 


AA 


Air A uT 


Av Aw 



Chis) 


Therefore, the force AF exerted on this surface element is given by 


AF ■=> f AA Ci'i^) 

where "p is the solar radiation pressure vector acting on the siirface 
element. Under very general conditions of surface reflection and 
re-radiation, it can be shown that this pressure vector has the form 

^ [C ^ i- C,)fi] 0-n) 


where S is the solar radiation fliuc per unit area normal to the flux, 
c is the velocity of light, and ^is the angle between the sun vector 
and the normal to the surface element, i.e., 

Cfi ^ ss n • Cf' 1^) 


The coefficients , Cg and may change with time due to aging of 
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of the surface material by some complex process. 

For the case of specular reflection and diffuse (Lambertian) 


reflection, the Cj^*s are given by 





(/• 2o) 

~ '3 

O' 2/) 


where «■ the fraction of incident radiation reflected specularly 
kg ” the fraction of incident radiation reflected diffusely 
by a Lambertian surface. 

It is to be noted that in equations (1.19) - (1*21), it is not 
implicitly assumed that the surface is radiating the entire energy 
incident on it, i.e., it is not necessary that we require the 
condition k^ + kg - 1 in order to obtain these equations. 

For the case of specul€ur reflection and non— Lambertian 
reflection and re -radiation a little consideration will reveal 

that the Cj^'s are given by 


W 

! 

C/'2U) 

c, = ^/«»- 

0'23j 

+ —O -0-r) 

1 

where / » the fraction of incident radiation reflected 

(specularly 

and otherwise) 



^ » the fraction of reflected radiation that is specular 
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- non-Lambertian coefficients for front and back 
surfaces respectively 

e^, e^ - emission values for front and back surfaces 
respectively. 

In passing, it may be noted that we have the relations 

/ y C/.3S) 

(I 

Moreover, it may be remarked that the form of equation (1.17) is 
valid for the more general non-Lambertian reflection and re-radiation 
which have a period of 7T in the azimuthal variable. In other words, 
Lambertian reflection means that the intensity I of the reflection 
is given by 

I :: 0'^7) 

Then, the case of hon-Lambertian reflection and re-radiation expressed 
by equation (1.2U) would correspond to an intensity which is indepen- 
dent of the azimuthal variable <j> and is of the form 

I ' I, fW 0.2S) 

where implicitly we exclude the case of Lambert’s law, i.e., 

^ c<n ^ 

The even more general case means that we can have reflection and 
re-radiation for which the intensity is of the form 

I - I, fee,'!’) O'io) 
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where 


f(e,f)=f(^>'t'+v) 0-^0 


Finally, to compute the total force F due to solar radiation, 
we obtain from equations (1.15) - (1.19) the following expressions 
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It is also trivial to compute the total torque M on the spacecraft 
by using the equations 

AM = ^ AF 0-1^) 

= 2 (axAF) 0-3i) 

but this will not be done here. 
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SECTION 3 - DISCUSSION 


It is obvious that the method just discussed does not 
encounter logic problems in determining self-shadovrlng . Moreover, 
because the points are first chosen on the projection plane, 

it circumvents the difficulties in book-keeping experienced in 
the other method of choosing vertices on the surface of the object. 
Furthermore, it does not require excessive core for storing the 
vertex data such as coordinates, area of surface element, nonnal 
vector, solar incidence angle, etc. This advantage becomes evident 
by evaluating the expressions in equations (1.32) - (1.3U) using 
three accumulators (one for each force component), not having to 
store the set of points . Finally, if greater accuracy is 

desired, it suffices only to choose smaller values Av^w*, multiply 

( A *A 

Av Aw j * perform 

computations only for the additional points newly introduced into 
the set ^v^,Wj|. This advantage cannot be realized in the other 
method of choosing vertices on the surface of the object. In that 
case, in going to a refined model with additional vertices, it is 
necessary to perform the entire computations starting from the 
beginning each time. 
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SECTION h " CONCLUSIOM 

From the foregoing discussion, it may be concluded that 
the present method has the following advantages: 

1, It does not experience logic problems in determining 
self-shadowing . 

2, It does not encounter the book-keeping problems arising 
in the case of choosing vertices on the svirface of the 
object. 

3, It does not require excessive core for storing vertex 
data. 

li. It can utilize previously obtained results in going to 
progressively more refined models. 
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